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ABSTRACT 

We introduce the Galactic Bulge Survey (GBS) and we provide the Chandra source list for the region 
that has been observed to date. Among the goals of the GBS are constraining the neutron star 
equation of state and the black hole mass distribution via the identification of eclipsing neutron star 
and black hole low-mass X-ray binaries. The latter goal will, in addition, be obtained by significantly 
enlarging the number of black hole systems for which a black hole mass can be derived. Further goals 
include constraining X-ray binary formation scenarios, in particular the common envelope phase and 
the occurrence of kicks, via source-type number counts and an investigation of the spatial distribution 
of X-ray binaries, respectively. The GBS targets two strips of 6° x 1° (12 square degrees in total), one 
above (1° < 6 < 2°) and one below (—2° < b < —1°) the Galactic plane in the direction of the Galactic 
Center at both X-ray and optical wavelengths. By avoiding the Galactic plane (—1° < b < 1°) we 
limit the influence of extinction on the X-ray and optical emission but still sample relatively large 
number densities of sources. The survey is designed such that a large fraction of the X-ray sources 
can be identified from their optical spectra. The X-ray survey, by design, covers a large area on 
the sky while the depth is shallow using 2 ks per Chandra pointing. In this way we maximize the 
predicted number ratio of (quiescent) low-mass X-ray binaries to Cataclysmic Variables. The survey 
is approximately homogeneous in depth to an 0.5-10 keV flux of 7.7x10"^'* erg cm~^ s~^. So far, 
we have covered about two-thirds (8.3 square degrees) of the projected survey area with Chandra 
providing over 1200 unique X-ray sources. We discuss the characteristics and the variability of the 
brightest of these sources. 

Subject headings: accretion; accretion disks — stars; binaries — X-rays: binaries 
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1.1. Multi-wavelength observations of X-ray sources 

X-ray observations are excellent probes of coronally 
active and accreting sources. Whereas studies of X- 
ray sources in our Galaxy have to a large extent fo- 
cussed on bright systems, investigations of fainter source 
class es have typically b een done in the Galactic Center 
(e.g. iMuno et al.|[2003[) . There, however, crowding and 
extinction make it more difficult to identify the correct 
optical and/or infr ared counterparts for large fractions 
of the sources (e.g. iMauerhan et al.ll200i^ . 

It is clear from previous surveys that multi-wavelength 
observations are vital for classifying faint X-ray sources 
since X-ray spectral information alone is rarely sufficient. 
Classification is important for various science goals as we 
outline below. 

The Galactic Bulge Survey (GBS) we present in this 
paper was designed to allow multi-wavelength observa- 
tions of the detected X-ray sources. The GBS consists 
of Chandra and optical imaging of two strips of 6° x 1°, 
one centered 1.5° above the Galactic plane and the other 
1.5° below the plane. We have chosen this area as the 
source density is still high, but by excluding |&| < 1° we 
avoid the regions that are most heavily affected by ex- 
tinction and source confusion (see Figure [T|). In Table [T] 
we provide a reference table for the acronyms used in this 
paper. 
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TABLE 1 

Reference table for the acronyms used in this paper. 



ACIS Advanced CCD Imaging Spectrometer 

AGN active galactic nucleus 

AM CVn star AM Canum Venaticorum star 

ASCA Advanced Satellite for Cosmology and Astrophysics 

BB black body 

BH black hole 

Brems Bremsstrahlung 

BSC bright source catalog 

CALDB calibration database 

ChaMPlane Chandra multi-wavelength plane survey 

CTIO Cerro Tololo Inter-American Observatory 

CV cataclysmic variable 

EoS equation of state 

CBS Galactic Bulge Survey 

GLIMPSE Galactic Legacy Infrared Mid-plane 

Survey Extraordinaire 

HMXB high-mass X-ray binary 

HR hardness ratio 

HRI high resolution imager 

ID identification 

IP intermediate polar 

LMXB low-mass X-ray binary 

NS neutron star 

PSPC position sensitive proportional counter 

RASS ROSAT all sky survey 

RS CVn star RS Canum Venaticorum star 

UCXB ultra-compact X-ray binary 

UKIDSS UKIRT Infrared deep sky survey 

UKIDSS/GPS UKIDSS Galactic plane survey 

VVV Vista Variables in the Via Lactea 

W UMa star W Ursae Majoris star 



The discovery of a large sample of X-ray sources will 
allow us to identify rare X-ray binaries such as quiescent 
eclipsing black hole (BH) and neutron star (NS) low- 
mass X-ray binaries (LMXBs) and ultra-compact X-ray 
binaries (UCXBs; LMXBs that have orbital periods <1 
hour). For dynamical mass measurements in LMXBs one 
needs to measure three parameters: the radial velocity 
amplitude of the companion star {K2), the ratio between 
the mass of the companion star and the NS or BH mass 
(q) and the inclination (?). A measurement of the rota- 
tional broadening of the stellar absorption lines (usini) 
combined with K2 gives this determination of q. The rea- 
son is that in Roche lobe overflow systems, like LMXBs, 
the companio n star is tidall y forced to co-rotate with the 
binary orbit (|Tassou]||1988l ). 

The system inclination can be determined through 
modelling of ellipsoidal variations caused by distortion of 
the c ompanion star using multi-color optical light cu rves 
Ce.g. lOrosz fc Hauschildtl l2Q00l ICantrefl et al.ll2010t ). In 
systems with favorable viewing angles, the (X-ray) 
eclipse duration can be use d to accurately determine 
the inclination (jHornd 119851 ) . Since the inclination is 
constrained by the geometry, mass measurements in 
eclipsing systems are independent of the modelling that 
lies behind inclinations derived from ellipsoidal varia- 
tions. Quiescent eclipsing systems are prime targets 
for optical mass measurements. Such mass measure- 
ments provide constraints on the NS equation of state 
(EoS) and on the dividing line between NS and BH 
systems (e.g. lOzel et aTl I2010D . Constraining the NS 
EoS remains one of the u ltimate goals for NS studies 
(|Lattimer fc PrakashI 1200 4! and references therein) and 
mass measurements of BH systems would improve our 



estimates of the stellar mass BH mass-distribution with 
implications for supernova and Gamma-ray burst mod- 
elling. Furthermore, the BH sample selected in quies- 
cence, as in the GBS, is not susceptible to potential X- 
ray outburst duty cycle based selection effects that oc- 
cur when selecting BHs after they have become X-ray 
bright and returned to quiescence, which is now com- 
mon practise. Finding eclipsing LMXBs with optical 
or near-infrared counterparts bright enough that opti- 
cal or near-infrared spectroscopy can provide accurate 
mass measurements is the principal goal of the GBS. 

1.3. Goals of the GBS: binary formation and evolution 

The second main goal of the GBS is to study the ori- 
gin and evolution of the population of X-ray binaries. 
By comparing the observed number of sources per source 
class (e.g. CVs and LMXBs) with those that have been 
predicted on the basis of population synthesis models, 
one can constrain binary evolution models, in partic- 
ular the nature of the common envelope phase. Most 
compact binary sources responsible for high energy phe- 
nomena went through one or two phases of common- 
envelope evolution, such as CVs, A M Canum Venatico- 
rum stars (AM CVns) and UCXBs (jTaam fc Sandquistl 
l2000f ). However, that phase is not yet well understood. 
Compact binaries have much lower orbital energy and 
angular momentum than the progenitor binary that con- 
tained giants (Paczynski 197(1). The binary semi-major 
axis is thought to shrink mainly during a phase of un- 
stable mass transfer and ejection, i.e., the spiral-in. If 
the outcome of this process is derived by assuming that 
the change in orbital energy is enough to eject the gi- 
ant's mantle, the predicted properties do not match the 
observations of double white dwarf binaries. These prop- 
erties can be matched with the assumption that the gi- 
ant's mantle is ejected, c arrying a fixed fraction of the to- 
tal angular momentum (jNelemans "etal|[2000l ). but this 
begs the question how the required energy is provided. 
A more complete theoretical description is required that 
takes into account both energy and angular momentum. 

From the properties of single radio pulsars it 
has been inferred that they receive a kick at 
birth, due to asymmet ries in the supernova explosion 
(|Lvne fc Lorimen Il994f ). This has important conse- 
quences for the evolution and propertie s of X-ray bina- 
ries (e.g. ' Brandt fc Podsiadlowskilll995l lKalogeral[T996[ ) 
so in principle X-ray binaries can be used to constrain 
the properties of the kick. In particular, the question is 
whether black holes also receive a kick at formation and, 
if so, what its properties are (e.g . I White fc van ParadiisI 
11996?, 'Jonker & Nelemaii^ [200l iWillems et al.l l20'05r 
iMillcr- Jones et al. 2009). 

Four main steps must be followed in binary p opulation 
synthesis modehng (jPostnov fc Yungelsonll200a ). First, a 
set of initial conditions must be chosen. Key parameters 
which must be set include the initial mass function, the 
initial distribution of mass ratios in binaries, the initial 
binary fraction, and the initial distribution of orbital pe- 
riods. Secondly, any evolution with time of these initial 
conditions such as due to changes in star formation rate 
should be taken into account. Next, a recipe for (binary) 
stellar evolution must be specified. The major uncer- 
tainties in this recipe include the aforementioned uncer- 
tainties in common envelope evolution; uncertainties in 
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the neutron star and black hole natal kick distributions; 
and uncertainties in mass and angular momentum loss. 
Finally, a recipe must be developed for converting accre- 
tion rates, typically calculated on timescales of centuries 
or longer, into instantaneous observables such as the X- 
ray luminosity. Key uncertainties in this step largely 
involve understanding the spectra of accreting sources at 
different accretion rates, and understanding how disk in- 
stabilities will affect the distribution of luminosities of 
those sources which undergo disk accretion. 

To make progress on these issues we envisage a two- 
pronged approach that includes detailed studies of indi- 
vidual systems on the one hand, and of the population on 
the other hand. Any viable evolutionary scheme must be 
able to reproduce the specific properties (such as compo- 
nent masses, orbital period, age, system velocity) of each 
observed individual system. Any viable evolution scheme 
must also reproduce the distributions of and correlations 
between these properties in the population of X-ray bina- 
ries. Observationally this requires the discovery of large 
homogeneous samples of CVs and (ultra-compact) X- 
ray binaries, and the detailed follow-up of a number of 
individual systems. 

The large number of sources will also allow us 
to investigate the spatial distribution of LMXBs and 
this provides input to the LMXB formation scenarios. 
iJonker fc Nelemani ()2004l ) found that there is a signif- 
icant excess of LMXBs at -10° < I < 0° w i th re - 
spect to 0° < ; < 10°. iWeidenspointner et all ()2008l ) 
found an asymmetric distribution of the 511 keV line 
emission using INTEGRAL data which they suggested 
was due to an asymmetric LMXB distribution (although 
see iBandvopadhvav et al.l [2009). However, any spa- 
tial non-uniformity in the Galactic birth distribution 
of LMXBs is thought to be washed out by the natal 
kick imparted on the NS in the LMXB forma tion mod- 
els pr esented by Ivan den Heuvell (|1983D and iKalogerg] 
(|1998D . On the other hand, the LMXB evolution- 
ary model, involving triple star evolution, discussed by 
lEggleton & Verbunt (1986) provides a channel for the 
formation of LMXBs without a large kick velocity. In 
general, except for formation scenarios involving either 
direct collapse to a BH or an accretion induced col- 
lapse, a Blaauw kick should be imparted regardless of 
the ty pe of compact object formed during the supernova 
event (lBlaauwlll961). Evidence for velocity kicks is found 
in th e z-distribution of LMXBs (vari Paradijs fc W hit^ 
119951: IJonker fc NelemansI 12004) and in the v e locity 
distribution of radio p ulsars (|Lvne fc Lorimeii 119941 : 
Hansen fc PhinnevI 119971 ) . An improved spatial distri- 
bution of LMXBs in the Bulge will help to test if such 
kicks happen commonly or rarely. 

1.4. Previous X-ray surveys of the bulge 

Previous X-ray surveys of the Galactic bulge provide 
some general expectations for our survey. The Advanced 
Satellite for Cosmology and Astrophysics (ASCA) ob- 
served the central region of our Galaxy with \l\ < 45° 
and |6| < 0.4° in the 0.7-10 keV band, albeit with a 
resolution of 3 ' which led to a lar ge number of unclas- 
sified sources (jSugizaki et al.ll200lD . The ASCA survey 
discovered 163 sources down to a flux of w 3 x 10^^^ erg 
cm~^ which as a group have properties in common 
with Cataclysmic Variables (CVs), high-mass X-ray bi- 



A, 
15 




Galactic longitude ( ) 

Fig. 1. — The large black plus white rimmed boxes indicate 
the two 6° X 1° GBS fields in Galactic coordinates. The thick 
white contour made-up of small circles each of a 14' diameter, in- 
dicate the area covered by our Chandra observations until 2010. 
The grey scale image and contours depict the total absorption in 
the Sloan i'-b and filter, A^i , estimated from the COBE dust maps 
llSchlegel et al... 1998 ). The contours are at A^l values of 2, 5 (both 
in dark grey), 10, 20, 50 and 100 (in white). 

naries (HMXBs), quiescen t LMXBs and Crab-lik e pul- 
sars. A recent paper by lAnderson et al.l (|2010D finds 
a number of massive stars within the ASCA Galactic 
plane sources. The XMM-Newton Galactic plane sur- 
vey covered 3 square degrees along the Galactic plane 
down to Fx (2-10 keV) 2 x 10"" erg cm'^ s-\ finding 
roughly one-third to be soft sources sugges tive of nearby 
coronally active stars (jHands et al.l I2004D . while two- 
thirds were hard, absorbed sources. The latter group 
is dominated by active galactic nuclei (AGN) but it also 
has a substantial Galactic component, especially between 
10"^^ > Fx > 10"^^ erg cm"^ s~^ Optical identifica- 
tions of 30 sources with Fx > 7 x lO"^'' erg cm^^ s^^ 
led to the classification of 16-18 coronally active st ars, 3 
CVs, and 2 LMXB candidates (jMotch et al.ll2010D . Fi- 
nally, spectroscopic follow up of medium-duration Chan- 
dra observations of i iumerous Galac t ic pla n e fields (the 
ChaMPlane survey, iGrindlav etldl l2005t iRogel et al.l 
[2006) finds roughly one-third to be soft sources, largely 
stars, and two-third to be hard sources, including many 
AGN and a few CVs. 

In this paper we describe the modelling that provides 
rough estimates for the number of sources that we ex- 
pect in the GBS (§2). The Chandra X-ray data, the 
X-ray source list, and analysis of the X-ray source prop- 
erties are described in §3. A brief outline of the multi- 
wavelength imaging data is given in §4, followed by a 
comparison of ROSAT sources in the GBS area with the 
C/ianrfra -discovered sources is §5. We end with a short 
summary (§6) and outlook (§7). 

2. SOURCE NUMBER ESTIMATES 

We used a simple method to estimate roughly the num- 
ber of X-ray sources in the Bulge area and based on this 
we devised an optimal observing strategy to best meet 
our goals. For each source class we used estimates of 
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their typical optical (i'-band) and X-ray brightness as 
well as local space densities or total number of objects 
in the Galaxy without following their formation or evo- 
lution. We distribute these sources in the Galaxy ac- 
cording to the model for the star formation history from 
[Nel emans ct al. (20(H), whic h is based on the Milky Way 
disk formation simulation of iBoissier fc PrantzosI ()1999f ) 
with an added Bulge component. We do not take the 
influence of kicks for neutron star systems into account. 
For each object we determine the interstellar absorption 
from the Schlcgel ct al. (1998) dust maps, assuming the 
dust is evenly distributed (i.e. homogeneous) between the 
Earth and the position where the line of sight leaves the 
dust layer for an assumed dust height above and below 
the Galactic plane of 120 pc. For the optical we use 
Av = i.lE{B - V) and A^, = Q.QAy. We further use 
Ak = 0.1 Ay. For the X-ray absorption we distinguish 
hard sources (with an assumed power law spectrum with 
photon index 2), soft sources (where the assumed spec- 
trum is a black body with temperature 0.25 keV), and 
sources with a Bremsstrahlung spectrum with a temper- 
ature of 2 keV. We link the absorption to the reddening 
via A^H = 0.179 x W^Ay (Prcdchl & Schmitt 1995). 

By comparing the expected number of quiescent 
LMXBs with the number of CVs as a function of X-ray 
flux (see Figure [2) we decided on rather shallow Chan- 
dra observations (« 2 ks per pointing) yielding a limiting 
flux of (1 — 3) X 10"^"* erg cm~^ s~^. This limit is depen- 
dent on the source spectrum including extinction. For 
the source number estimates we take a general flux limit 
of SxlO"^** erg cm~^ s~^. Deeper observations would 
mostly increase the number of detected CVs. This would 
make it more difficult to classify the LMXBs among the 
sample of X-ray sources. As mentioned above, we se- 
lected two strips of 6° x 1°, one centered at 1.5° above 
and one centered at 1.5° below the Galactic Center (see 
Figured]). The size and location are chosen in order to 
detect enough sources in each class to achieve our science 
goals. The optical extinction is still substantial in this re- 
gion but relatively low compared to the Galactic Center 
and Galactic plane. This will allow us to obtain optical 
spectroscopic follow-up to a large fraction of the X-ray 
sources, which is crucial for achieving our science goals. 
Going to a higher latitude would facilitate identification 
as the crowding would be less severe, but the source den- 
sity drops off quickly, which would require a larger sur- 
face area to detect a sufficient number of sources. The 
chosen strips are thus a compromise between survey area 
size and location. 

The number of sources expected in the GBS area and 
the assumptions made in the population estimates are 
given in Table [2J The values for the assumed X-ray 
luminosity, i'-band absolute magnitude and optical ex- 
tinction Ay in Table [2] are averages. In the modelling 
we assumed a distribution of X-ray luminosities using 
a Gaussian distribution, with cr = 0.5 and the average 
normalized to 1. For the i'-band absolute magnitude we 
apply a Gaussian smoothing with a—\ magnitude. Fi- 
nally, we use a Gaussian smoothing in the assumed Ay 
with fj — Q.2 magnitude. 

We expect to discover Roche lobe overffow NS and 
BH X-ray binaries, both with main sequence donors 
(LMXB) as well as UCXBs. Most LMXBs wfll be in 
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Fig. 2. — The predicted number of non-magnetic CVs (thick blue 
lines), intermediate polars (thin red lines) and quiescent LMXBs 
(medium thickness black lines) in the GBS area as a function 
of source X— ray flux. The approximate flux limit of the GBS is 
(1 — 3) X 10~^ erg cm~-^ s~^, with the variations largely due to 
differences in source spectra. We expect that deeper Chandra ob- 
servations would increase the number of CVs much more strongly 
than they would increase the number of quiescent LMXBs. 



quiescence. The numbers for the LMXBs are based on 
a very simple estimate of «140 p ersistent systems w ith 
Lx > lO^^ergs-i in the Galaxy (jCrimm et al.ll2002D . a 
number ratio of quiescent to persistently bright/active 
LMXBs of ~ 70, and the prediction that the number of 
UCXBs r riay be equal to the numb er of non-degenerate 
LMXBs (|Belczvnski fc TaamI [200l . To estimate the 
number ratio of quiescent to active LMXBs we use re- 
sults of surveys of Galactic globular clusters as input. 
There, roughly 10 times as many quiescent LMXBs are 
identified (by their soft, blackbo dy-like spectra) as ac- 
tive LMXBs (|Heinke et all l2003| ). Since fainter quies- 
cent LMXBs tend to be dominated by their hard, power- 
law component rather than the blackbod y component 
(e.g. iJonker et all 120041 and iJonkeii I2008D . it has been 
suggested that half of all quiescent LMXBs are missed in 
these globular cluster surveys (Heinke et al. 2005) . How- 
ever, it is uncertain whether globular cluster LMXBs are 
similar in their duty cycles to Galactic LMXBs, as they 
are formed in different ways. 

An alternative way to estimate the number ratio be- 
tween active and quiescent LMXBs comes from the total 
numb er of LMXBs and the estimat ed X-ray binary life- 
time. iPortegies Zwart et al.l (|1997| ) compute that there 
are 10~^ NS binaries formed per year for most reasonable 
values of the common envelope parameter. Assuming 
a typical lifetime as an X-ray binary of approximately 
1 Gyr, then t here are « 10^ syst ems in the Galaxy. 
More recently, iKiel fc HurlevI ()2006( ) find a similar num- 
ber from their evolutionary population synthesis calcula- 
tions. Taking that there are «140 active X-ray binaries 
at any given time, one in seventy are active at any given 
time. We take this 70 as our assumed ratio between ac- 
tive and quiescent LMXBs in our number estimates. 

In an X-ray selected sample such as this, another ma- 
jor population of sources that we will pick-up is CVs 
and of these we will in particular detect magnetic In- 
termediate Polar systems (IPs). Their number esti- 
mate is based on the observed IP and non-magnetic 
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TABLE 2 

Estimated numbers of the various source categories. The columns describe: the source type (I), the assumed mean X-ray 

LUMINOSITY IN ERG S~l (II), THE X-RAY COLOR (III), THE ASSUMED MEAN i'-BAND ABSOLUTE MAGNITUDE (IV), ASSUMED i' — K 
INTRINSIC COLOR (V), ASSUMED TOTAL NUMBER OF SOURCES IN THE GALAXY (INTEGER VALUES) OR THE SPACE DENSITIES (p PC"^) (VI), 

NUMBER OF SOURCES THAT SHOULD BE DETECTABLE BOTH IN THE GBS X-RAY AS WELL AS THE OPTICAL SURVEY (VII), NUMBER OF 
SOURCES DETECTED IN THE GBS X-RAY AS WELL AS IN A K-BAND IMAGING SURVEY OF LIMITING MAGNITUDE OF = 18 (VIII), NUMBER 
OF SOURCES DETECTED IN X-RAYS IN THE GBS AREA (IX), NUMBER OF SOURCES FALLING IN THE GBS AREA (X). QUIESCENT SOURCES ARE 
DENOTED WITH 'Q'. 'BB' STANDS FOR A 0.25 KEV 'BLACK BODY'. 'BrEMS' STANDS FOR BrEMSSTRAHLUNG. 



(I) 


(II) 


(III) 


(IV) 


(V) 


(VI) 


(VII) 


(VIII) 


(IX) 


(X) 


LMXB 


IQ35 


Hard 








140 


6 


7 


7 


7 


qLMXB 


Iq33 


BB 


5 


2 


10000 


120 


86 


221 


532 


UCXB 


1034 


Hard 


4 





1000 


32 


3 


56 


58 


qUCXB 


1032 


Hard 


10 





10000 


1 





8 


605 


CV (non mag.) 


1031 


Brems 


7.5 





2 X 10-5 


62 


61 


62 


1.4 X 10^ 


CV (IP) 




Brems 


8.5 





1.5 X 10-** 


152 


5 


525 


7.7 X 104 


RS CVn 


1031 


Hard 


2.5 


1 


1 X 10-4 


596 


596 


596 


1.3 X 10'^' 


W UMa 


5 X 1030 


Hard 


4.5 


2 


7.5 X 10-5 


160 


160 


160 


2.3x10''^ 


Be X-ray binaries 


1034 


Hard 








500 


9 


9 


10 


10 


Total 












1142 




1648 





CV sp ace densit ies of ~ 1-5 x 10 ^ and 2 x 10 ^ pc 
(iPatter son 198 ^Hertz et al.lll990L iPretorius et"aIll2007L 
iRogel et al.i i2008D . T he discovery of a lar ge number of 
IPs with INTEGRAL (jScaringi et al.ll2010D may be indi- 
cating that the fraction of IPs among the CVs is substan- 
tially larger than was previously thought. In addition, 
there will be a large number of active (binary) stars such 
as RS CVn sta rs. The expected nu mljer of background 
AGN is 25 (see lEbisawa et al.ll2005l ). We might also find 
a few ultra-compact double white dwarf systems (AM 
CVns). 

The expected number of W UMa-like sources is cal- 
culated using the space density of 7.5x10-^ pc"'^ with 
Mi' fti 4.5 and an av erage X-ray luminosity of 5 x lO'^'^erg 
s-i H ucinskil[l998[ ). Finally, the expected number of Be 
- X-ray binaries is estimated using an estimated total 
number of systems of 500 in our Galaxy, with an Mj' = 
and a hard power law spectrum with an average luminos- 
ity of Lx = 1 X lO^^erg s'^ We expect 10 in the GBS 
area out of which 9 will be detectable in the i'-band as 
well as the isT-band. 

These numbers are uncertain and calibrating the num- 
ber densities of sources using the GBS is one of the prin- 
cipal objectives. We will compare our identifications with 
the predicted numbers of binaries in each category and 
thus place strong constraints on important uncertainties 
in the theory of binary evolution, such as the common- 
envelope phase and the existence and magnitude of kicks 
imparted on NSs and BHs. 

Of the predicted LMXBs we assume that 10 per cent 
will have BH accretors. This number is highly uncertain 
but 10 per cent is roughly the number rati o between NS 
and B H systems observed in HMXBs (see I Grimm et al.l 
120031 and references therein). Some authors have claimed 
that the rat io between NS and BH systems should be 
closer to 1 (|Romanil Il992f) , whereas others indeed find 
values close r to 10 per cent (jPortegies Zwart et al.lll997l : 
iKiel fc Hur lcv 200q). Assuming 10 per cent would lead 
to 25 BH LMXBs in the GBS which will approxi- 
mately double the population of known Galactic BHs for 
which a dynamical mass measurement is possible (cur- 
rently ~20, [McQintockiLRsmillajo^ Eclipsing 
BH LMXB systems should exist, but have not yet been 



found. It has been proposed that they are too weak to be 
detected by current X-ray all sky monitors because they 
are ob scured behind the accretion d isk rim during out- 
burst (|Naravan fc McClintocgi2005D . If so, they should 
turn up in our GBS in quiescence, when the disk is ex- 
pected to be much thinner. The number of eclipsing 
sources depends on the distribution of the mass ra tio be- 
tween the accretor and the donor star (see e.g. iHornd 
11985( 1. For mass ratios q ~ 0.3 approximately 20-25 per 
cent of the 200 new quiescent LMXBs we expect to dis- 
cover should be eclipsing of which a handful could be 
BHs. 

Overall, this modelling can be summarized in a figure 
with X-ray flux as a function of optical «'-band magni- 
tude for the main components of the expected sources 
(Figure [3]). The large majority of bright optical sources 
with i'< 16 associated with an X-ray source in the GBS 
area is expected to be an RS CVn. Similarly, the major- 
ity of X-ray bright sources (Fx > 1 x 10^^^ erg cm^^ s^^ 
with a optical counterpart with 18 < i' < 22 is expected 
to be an active LMXB or UCXB. Finally, most quiescent 
LMXBs and UCXBs will have faint optical counterparts 
[i' > 21). 

Clearly, other source types are expected, for in- 
sta nce single X-ray activ e nearby G, K or M stars 
(cf. lSchmitt fc Liefkd[200l and Algol sources. However, 
the expected number of sources is small for the Algol 
and nearby single stars, so we do not provide detailed 
estimates for these classes. 

In order to roughly assess the accuracy of the afore- 
mentioned modelling, we compare the number of sources 
detected in the single, 1 Ms deep Chandra observation of 
a small part of 20.6 square arcminute s of the GBS area 
presented by iRevnivtsev et al] ()2OO90 with the number 
of sources we predict if one were to observe the whole 
GBS area to that depth of 1 x 10~^^ erg cm"^ s^^. We 
expect to find 3 million sources based on our source num- 
ber estimates, whereas scaling the observed 20.6 squared 
arcminutes to 12 square degrees we would expect to find 
about 1 million sources. The main difference between 
these two numbers probably stems from uncertainties 
in scaling from the 20.6 square arcminute area to a 12 
square degree field such as those caused by the differ- 
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Fig. 3. — Predicted source X-ray flux, Fx, as a function of i'-band magnitude for the sources predicted to be most numerous in our 
GBS survey. The vertical dashed line indicates the approximate i'-band magnitude limit of our optical imaging. 



ence in extinction between th e GBS area with respec t 
to that in the area studied bv iRevnivtsev et all (|2009( ). 
as weh as of course uncertainties in our modelhng. Ah 
in ah, the agreement within a factor of 3 between the 
number of sources we predict and what was found is not 
unsatisfactory. 

3. X-RAY OBSERVATIONS 

We have obta ined observations with the Chandra X- 
ray observatory (Weiss kopf et aLll2002| ) covering, to date, 
about two-thirds of the total area of twelve square de- 
grees that we envisage for the GBS. 

In Figure [TJ the wiggles indicate the composite outline 
of each circular field of view of 14' diameter of the indi- 
vidual Chandra observations obtained to date covering 
the GBS area. The Chandra observations have been per- 
formed using the 10 to 13 CCDs of the Advanced CCD 
Imaging Spectrometer (ACIS) detector (jCarmird 119971 : 
ACIS-I). The observation identification (ID) numbers for 
the data presented here are 8643-8774 and 9977-10024. 
We reprocessed and analyzed the data using the CIAO 
4.2 software developed by the Chandra X-ray Center 
and employing CALDB version 4.3. The data telemetry 
mode was set to very faint for all observations except that 
with ID 8687 since there was a bright source present in 
archival ROSAT observations. For that observation we 
used the standard faint mode. 

We also reprocessed and reanalyzed the data that 



was obtained as pa r t of the Bulge Latitude Survey 
(jCrindlav et all 120051 : iHong et all 120091 observation IDs 
7160-7162, 7166-7168, 8199-8204, 9562-9564) that falls 
in the area we target in the GBS using exactly the same 
reduction and analysis as for our GBS observations. Be- 
cause the original Bulge latitude survey Chandra obser- 
vations have an exposure time of approximately 15 ks, 
we selected 2 ks segments to allow for a comparison with 
our GBS observations. We selected data stretches of 2 
ks length taking the start time of the observation from 
the header of the data, plus 100 seconds as the starting 
point of the 2 ks stretches. These data were also obtained 
using the very faint mode. 

The very faint mode provides 5x5 pixel information 
per X-ray event. This allows for a better screening of 
events caused by cosmic rays. In our analysis we selected 
events only if their energy falls in the 0.3-8 keV range. 

We used wavdetect to search for X-ray sources in 
each of the observations using data covering the full 0.3- 
8, the 0.3-2.5 and the 2.5-8 keV energy band separately. 
We set the sigthresh in wavdetect to 1x10"^, which 
implies that for a background count rate constant over 
the ACIS-I CCDs there would be 0.1 spurious source 
detection per observation as about 1 x 10^ pixels are 
searched per observation. However, as we explain below, 
we applied additional selection criteria. This lowers the 
number of spurious sources. 

The resulting SIGNI column in the output hst of de- 
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tected sources provided by wavdetect is an estimate of 
the photon flux significance, not of the detection signif- 
icance. Instead, we retained all sources for which Pois- 
son statistics indicates that the probability of obtaining 
the number of detected source counts by chance given 
the expectation for the local background count rate is 
lower than 1x10^®. This would be equivalent to a > 5cr 
source detection in Gaussian statistics. Next, we deleted 
all sources for which wavdetect was not able to provide 
an estimate of the uncertainty on the right ascension [a] 
or on declination [S] as this indicates often that all counts 
fell in 1 pixel which could well be due to faint afterglows 
events caused by cosmic ray hits. In additi on, we impose 
a 3 co unt minimum for source detection as lMurrav et al.l 
(|2005f) simulated that in their XBootes survey with 5 ks 
ACIS-I exposures, 14 per cent of the 2 count sources 
were spurious (note that this percentage will probably 
be lower for our GBS exposures of 2 ks). 

Since our Chandra observations were designed to over- 
lap near the edges, we searched for multiple detections 
of the same source either in one of the energy sub-bands 
or in the full energy band. We consider sources with 
positions falling within 3" of each other likely multiple 
detections of the same source. Finally, we inspected the 
source list and found two spurious sources caused by the 



bright read-out trail and piled-up core of source # 1 (see 
Table [S]). 

In total we detected 1234 distinct sources in the area 
indicated with circles in Figure [TJ including sources de- 
tected in the Bulge Latitude Survey area that fall within 
our GBS area. The source list is given in Table [3] and 
provides information on a, S, the error on a and 6, total 
number of counts detected, the observation ID of the ob- 
servation resulting in the detection and the off-axis angle 
at which the source is detected. The error on a and S 
are the error provided by wavdetect, it does not take 
into account the typical Chandra bore-sight uncertainty 
of 0.6" (90 per cent confidence). 

We provide individual Chandra source names, how- 
ever, for briefness we use the source number in Ta- 
ble |3] to indicate which source we discuss in this pa- 
per. For the error crjv on the detected number of counts 
N, iGrimm et al.| (|2^ give aN = I + VN + 0.75 af- 
ter iGehrcls (198a). To allow for an rough, easy calcula- 
tion of the source flux based on the detected number of 
source counts we give the conversion factor for a source 
spectrum of a power law with photon index of 2 ab- 
sorbed by TVh = 1 X 10^2 cm-2: 7.76 x 10 " 
s^^ photon^^. 
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TABLE 3 PLACEHOLDER, FIRST TEN ENTRIES ONLY! See 
|htt p: / /www. sron.nl/~pcterj/gbs or the electronic version in ApJS for the full 
source list. The UBS X-ray source list providing the GBS source name, the 
source number as used in this paper, a, S, the error on a and S, total number 

of counts detected, the observation ID of the observation resulting in the 
detection, the off-axis angle at which the source is detected and the hardness 

ratio (HR) for sources detected with more than 20 counts. The HR is 
calculated for the detection where the off-axis angle was smallest if the source 
was detected multiple times. 
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CXOGBS 


J175024, 


.4- 


■290216 


1 


267.60182588 


-29.037885415 


0.034 


0.011 


3391 


8709 


8.99 


1 




CXOGBS 


J173728, 


.3- 


■290802 


2 


264.36831021 


-29.133892744 


0.058 


0.084 


2191 


8691 


7.63 


2 




CXOGBS 


J174042, 


.8- 


■281808 


3 


265.17839123 


-28.302224391 


0.048 


0.033 


1850 


8687 


5.54 


4 




CXOGBS 


J173931, 


.2- 


■290952 


4 


264.88008702 


-29.164675805 


0.104 


0.072 


238 


8679 


5.63 


1 


-0.78 ±0.09 


CXOGBS 


J174009, 


.1- 


■284725 


5 


265.03805780 


-28.790455138 


0.055 


0.054 


157 


8677 


3.41 


1 


0.50 ±0.10 


CXOGBS 


J174445, 


.7- 


■271344 


6 


266.19074704 


-27.229022849 


0.065 


0.047 


153 


8647 


3.20 


1 


-0.25 ±0.09 


CXOGBS 


J173826, 


.1- 


■290149 


7 


264.60910401 


-29.030389506 


0.160 


0.155 


150 


8690 


6.65 


2 


-0.89 ±0.12 


CXOGBS 


J173508, 


.2- 


■292957 


8 


263.78448781 


-29.499426909 


0.119 


0.073 


138 


9997 


4.78 


1 


-0.39 ±0.10 


CXOGBS 


J173508, 


.3- 


■292328 


9 


263.78498489 


-29.391224268 


0.141 


0.135 


134 


9996 


5.95 


1 


-0.64 ±0.12 


CXOGBS 


J173629, 


.0- 


■291028 


10 


264.12099622 


-29.174663351 


0.322 


0.180 


122 


9995 


8.34 


3 


-0.76 ±0.16 



3.1. Sources detected multiple times 

In the trade-off between a homogeneous survey depth 
and total survey exposure time we used 7' as an effec- 
tive radius of the Chandra field of view in designing the 
survey. The size and shape of the point spread function 
for off-axis angles larger than 7' has degraded such that 
many optical and/or near-infrared stars will fall inside 
the X-ray error circle, even for sources detected with 
more than 10 counts , making multi-w avelength follow- 
up more difficult (cf. iHong et al.ll2005l ). As noted above, 
this yields the possibility of multiple detection of sources 
discovered in the overlap regions. Indeed, 105 sources 
are detected more than once, where we have taken source 
positions within 3" of each other as multiple detections 



of the same source. Out of these 105 sources, 95 sources 
are detected two times, 9 sources are detected three times 
and 1 source is detected 4 times. The properties that we 
list in Table [3] for these sources are those of the detection 
that gave rise to the largest number of X-ray counts. In 
Table [3] we also list the number of times that sources are 
detected. In addition, we found that although Chandra 
source # 33 and # 230 are 4.8" apart, they are con- 
ceivably two detections of the same source. Both detec- 
tions are rather far off-axis (5.2' and 8.8', respectively) 
explaining the relatively large uncertainty in both mea- 
surements of the position. 

3.2. Detection probability vs off-axis angle 

We investigate the relationship between the number of 
sources detected as a function of off-axis angle and the 
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Fig. 4. — The number of detected sources normalized to area on 
the sky as a function of off-axis angle. The black solid line in- 
dicates all sources, the red dashed line indicates sources detected 
with only 3 X-ray counts and the blue dot-dash line are sources 
detected with 10 counts or more. From the black and red lines it 
is clear that the survey sensitivity to low count sources is higher 
when nearer to on-axis, as expected. The fact that the blue line is 
approximately flat up to 7' indicates that the GBS survey sensitiv- 
ity is approximate constant for sources of 10 counts or more. For a 
power law spectrum with photon index of 2 and an Njj = 1 X 10^^ 
cm~^ this corresponds to a flux of 7.7x10"^* erg cm~^ s"'^. 

number of source counts. In this way we can determine 
the approximate number of source counts where the sur- 
vey sensitivity for off-axis angles less than 7' is constant. 
We find (Figure |4]) that when normalized to the surface 
area, the number of detected sources drops as a function 
of off-axis angle for angles larger than «2'. The prob- 
able reason for the slightly lower number of sources de- 
tected per unit surface area inside the 2' radius is that 
the chip gaps between the 4 ACIS-I CCDs subtend a 
larger fraction of the total solid angle within 2' than 
they do in the outer regions. That area has a lower ef- 
fective exposure time as can be seen in Figure [51 For 
source counts larger than 10 the normalized number of 
detected sources is approximate constant for off-axis an- 
gles less than 7' (blue dot-dash line in Figure S]). Thus 
the GBS coverage is approximately homogeneous down 
to source fluxes of w 7.7 x 10~^^ erg cm~^ s~^, whereas 
the sensitivity is decreasing with off-axis angle for fainter 
sources. 



3.3. X-ray spectral information 

We extract source counts using circular source extrac- 
tion regions of 10". Background extraction regions are 
annulli with inner and outer radii of 15" and 30", respec- 
tively. We plot the 89 sources for which we detected more 
than 20 counts in a hardness - intensity diagram (Fig- 
ure H]) . To mitigate the effects that small differences in 
exposure time across our survey can have, we use count 
rates as a measure of intensity. We define the hardness 
ratio as the ratio between the count rate in the 2.5-8 
keV minus that in the 0.3-2.5 keV band to the count 
rate in the full 0.3-8 keV energy band (after iKim et all 
[2004) . We derived the hardness using XSPEC version 




Fig. 5. — The exposure map for the observation with Obs ID 
9999 for an assumed mono-chromatic source emitting 1 keV pho- 
tons, showing the lower exposure for the gaps between the ACTS I 
CCDs. Lighter areas have a higher exposure on the sky (in terms of 
cm^ Due to the satellite dithering some exposure is achieved 

even for the gaps. The vertical strips are caused by the presence 
of columns of pixels of reduced efficiency. The grey area is approx- 
imately 8' X 8' in size. 



12.6 ([Arnaudl 119961 ) by determining the count rates in 
the soft and hard band taking the response and ancillary 
response file for each of the sources. The three brightest 
sources suffer from photon pile-up, where more than one 
photon is registered during the CCD integration time 
(see e.g. IDavisI [2OOII) . Photon pile- up will artificially 
harden those 3 sources in a hardness-intensity diagram, 
therefore, we have not plotted these 3 sources in Fig- 
ure [HI nor in Table [3| For the other 86 detected sources 
photon pile-up is unimportant (less than a few percent). 
Naively, one would expect most hard sources to be more 
distant and more reddened than the soft sources, as the 
intrinsic spectral shape of the most numerous classes of 
sources we expect to find does not differ much. 

The most interesting aspect from Figure [H] is perhaps 
the paucity of bright hard sources. Conversely, most 
bright sources are probably nearby with low extinction. 
As foreseen, the spectral information is insufficient for 
source classification for the majority of the total number 
of detected sources, therefore, classification will have to 
come from multi- wavelength observations. Finally, there 
seems to be a dichotomy in the hardness with one peak 
centered on a hardness of 0.4 and another centered on 
-0.5 with a clear paucity of sourc es with hardne s s 0. A 
similar dichotomy was reported in [Warwick et al.l (j2010D . 

A conceivable explanation for the nature of this di- 
chotomy is that the soft sources are dominated by nearby 
sources, such as RS CVns. The harder sources are further 
away, and are either intrinsically hard spectrum sources 
or they appear hard due to the effects of extinction. The 
latter increases strongly towards the Gal actic cen ter for 
distances between « 3 — 5 kpc (Marsh all et al.l [2006) . 
The strong increase in extinction over a small distance 
interval would, in this scenario, be responsible for the 
lack of sources near a hardness of 0. 

4. MULTI- WAVELENGTH OBSERVATIONS 
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Fig. 6. — The hardness — intensity diagram for the 86 sources 
for which more than 20 counts were detected in the GBS survey 
(we exclude the three brightest sources from the plot since their 
hardness ratio is strongly affected by effects of photon pile-up). 
To mitigate effects of small differences in exposure times we used 
count rates as a measure of intensity. The hardness is defined as 
the ratio between the count rate in the 2.5-8 keV minus that in 
the 0.3-2.5 keV band to the count rate in the full 0.3-8 keV energy 
band. Ifard sources fall in the top half and soft sources in the 
bottom half of this figure. The red line with large bullets shows 
the influence of the extinction Njj on a power law spectrum with 
index 2 for a source count rate of 0.05 counts s~^ and Njj of 0.01 
xlO^^ cm~^. The Njj is increasing from bottom right to top left 
from (0.01, 0.1, 1, 3, 10)xl022 cm-2. 

Optical imaging observations covering the full GBS 
area have been obtained using the prime focus MO- 
SAIC II instrument mounted on the 4 m Victor 
M. Blanco telescope at the Cerro Tololo Inter- American 
(CTIO) observatory. We employed the Sloan ; i, and 
Ha filters with exposures of 120, 180 and 480 seconds, 
respectively. 

For the optical astro metry we used the UCAC2 catalog 
(jZacharias et al.ll2004[ ). The UCAC2 is made up by stars 
brighter than 16^^ magnitude that would saturate in the 
120 s-long --band observations. Therefore, to facilitate 
accurate astrometric calibration of the optical observa- 
tions we have also obtained short, 10 s long, --band ex- 
posures for each pointing. A detailed discussion of the 
optical observations and the findings will be presented in 
a forthcoming paper (Bassa et al. in preparation). 

Using the same telescope and instrument (Blanco and 
MOSAIC II), we re-imaged the fields of the initial opti- 
cal observations at random intervals over 9 nights in the 
Sloan --band. The goal of these observations is to search 
for variability induced by binary motion. A detailed dis- 
cussion of these optical variability observations and the 
findings will be presented in a forthcoming paper (Hynes 
et al. in preparation). 

In addition to the optical observations, the survey 
field has also partially been observed in the near- 
infrared as part of the UKIRT Infrared deep sky sur- 
vey (UKIDSS) G alactic plane survey (UKIDSS/GPS; 
iLucas et al.l 120081 ) . Furthermore, the GBS area is fully 
covered as part of the Vista Variables in the Via Lactea 



(VVV) survey (|Minniti et al.l l2010t Greiss, Steeghs et 
al. in prep). These near-infrared observations will also 
help us distinguishing between the different expected 
source types. Finally, the Spitzer Galactic Legacy 
Infrared Mid-plane Survey Extraordinaire (GLIMPSE; 
[Churchwcll et al .. 2009f ) survey has covered the GBS area, 
providing mid-infrared imaging of the survey area. 

We have also started with the optical spectroscopic 
identification of the detected X-ray sources using a suite 
of optical telescopes such as the Blanco telescope, the 
New Technology Telescope, the Magellan telescopes, the 
Gemini-South telescope and the Very Large Telescope. 
Results of these observations will also be presented in 
forthcoming papers. 

5. COMPARISON WITH ROSAT ALL SKY SURVEY 
SOURCES 

In order to investigate whether bright sources in our 
source l ist are detec ted by the ROSAT All Sky Survey 
(RASS: rVoges et aLiiiggg ) and, conversely, to investigate 
if bright RASS sources are still detected in our Chan- 
dra observations, we searched for ROSAT sources within 
30" of our Chandra positions. 

5.1. Persistent sources among Chandra source #1-10 

The brightest Chandra source, source ^1 is a 
transient (see below). The four ROSAT Bright 
Source Catalog (BSC) sources, IRXS J173728.0- 
290759, IRXS J174043.1-281806, IRXS J173933.4- 
291001, IRXS J173826.7-290140, probably correspond to 
Chandra source number 2, 3, 4, and 7 in Table [31 as the 
nominal offset between the Chandra and ROSAT posi- 
tion is, 5.76, 4.26, 29.7, and 11.6", respectively. The 
ROSAT positional (la) uncertainty of these four sources 
is 9, 8, 13, and 16", respectively. Source #2 of the 
Chandra GBS survey (IRXS J173728.0-290759) is a well- 
known Seyfert 1 galaxy Source #3 (IRXS J174043.1- 
281806) is the persistent LMXB and UCXB-candidate 
SLX 1737-282. Source #4 is associated with HD 316072; 
the angular distance between the Chandra position and 
HD 316072 is 0.46". HD 316072 is a bright, V=9.92, 
G9III star. This source is conceivably a bright ac- 
tive star or a long-period LMXB. Taking the observed 
and absolute y-band magnitude and assuming Ay = 
yields a maximum source distance of 675 pc. Mod- 
elling the Chandra X-ray spectrum consisting of 238 
photons, we find that the source spectrum is soft. It 
can be described by an absorbed black body spectrum 
with kT=0.22 keV and Nh = 0.4 x 10^^ cm^^. This pro- 
vides a good fit at an absorbed/unabsorbed 0.5-10 keV 
flux of 5.6xl0~^'^/1.2xl0~^^ erg cm~^ s~-^, respectively. 
This converts to a 0.5-10 keV luminosity upper limit of 
6.5xl0'^^erg s~^. Such a soft spectrum would be con- 
sistent both with a quiescent NS LMXB as well as a 
white dwarf accretor in a CV. Optical time series spec- 
troscopy should allow us to distinguish between an active 
(binary) star or an X-ray binary scenario. Source #7 is 
0.34" away from a bright optic al sourc e identified as a 
pre-main sequence star bv iTorres et al.l (p006i ). 

Similarly, sources #5, 6, 8-10 are bright enough that 
they should have been detected by the RASS BSC. 
Source #5 is 2.78" away from AX J1740. 1-2847, which 
is identifi ed as a low -luminosity high-mass X-ray binary 
pulsar (|Kaur et al.lfeOlO.) . Source #6 is a known Be X- 
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ray binary which is hsted in the RASS faint source cat- 
alog as IRXS J174444. 7-271326, indicating some moder- 
ate X-ray variabihty. Source #8 is nominahy 29.8" away 
from AX J1735. 1-2930, which is currently unclassified. 
Source #9 is 0.33" away from HD315997, which is an 
A5 star in a known eclipsing binary in a 2.8723 d orbital 
period (jNesterov et al.lll995l: lOtero et al.ll2006[ ). The V- 
band magnitude varies between 11. 21-11.40 with a sec - 
ondary eclipse minimum of 11.24 (lOtero et all 120061) . 
The small size of a neutron star and black hole implies 
that an accretion disk should be present if one wants 
to explain the amplitude and duration of the secondary 
eclipse in an X-ray binary context. Taking the observed 
and absolute F-band magnitude and assuming Ay = 
yields a maximum source distance of 695 pc for the A5 
star. This implies a maximum source 0.5-10 keV lumi- 
nosity of 6 X lO^^erg s~^. It is also possible that the 
A5 star is in orbit with a late type star that is not seen 
in the optical spectrum but which is responsible for the 
X-ray emission. Again, orbital phase resolved spectro- 
scopic observations are necessary to reveal the nature of 
this object. 

Source #10 is 0.47" away from HD 315992 
(F8, V=9.98) and nominally 26.38" away from 
IRXS J173628. 8-291055 which is listed in the RASS faint 
source catalog. In the ASCA Galactic center survey the 
source is f ound to have a 0.7- 10 keV flux of lxlO~^^ erg 
cm~^ s~^ (|Sakano et al.l[2002l ). As for several of the other 
sources discussed above, multiple optical spectra of this 
source will allow us to investigate if they are in a binary 
and, if so, what the nature of the second star is. 

To summarize, the Chandra point sources #2-10 cor- 
respond to: the Seyfert 1 galaxy IRXS J173728.0- 
290759, SLX 1737-282, HD 316072, AX J1740.1- 
2847, IRXS J174444.7-271326, IRXS J173826.7-290140, 
AX J1735. 1-2930, HD 315997 and HD 315992. 

5.1.1. Chandra light curves 

We inspect the Chandra light curves of source #1-10. 
Source #1, 3, and 5 show evidence for flare-like variabil- 
ity. Fitting the light curve with a constant gives a 
value of 359, 57, and 43 for 20, 17, and 20 degrees of 
freedom, respectively. 

There is marginal evidence in the light curve of source 
#4 for the presence of an eclipse-like feature (see Fig- 
ure [7]), although, fltting the light curve with a constant 
gives a value of 27.5 for 20 degrees of freedom. 

The light curves of source #2, 6, 7, 8, 9, and 10 are 
consistent with being constant with values of 24.3, 
27.2, 22.7, 27.4, 12.4 and 20.5 for 20 degrees of freedom, 
respectively (except source #9 for which there are 19 
degrees of freedom). 

5.2. Transient sources 

The only other two BSC RASS sources in the GBS 
area, IRXS J175113.2-293842 and IRXS J174220.8- 
273736, are not detected in our GBS survey. Their nom- 
inal ROSAT positional errors are 26" and 12", respec- 
tively. This means that, especially for IRXS J175113.2- 
293842, it is not inconceivable that the source has a 
Chandra counterpart further away than 30". However, 
even in a search radius of 1' there is no Chandra source 
possibly related to IRXS J175113.2-293842. Whereas 
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Fig. 7. — The light curve with a 200 s bin size of source #4 from 
Table [3] There is marginal evidence for a dip/eclipse like feature 
near T=1500 s, although the value of a fit of a constant is 27.5 
for 20 degrees of freedom. 

the band passes of ROSAT (0.1-2.5 keV) and Chan- 
dra (0.3-8 keV) are different, Chandra is much more 
sensitive and both sources should have been exception- 
ally soft for it not to be detected in Chandra; this is 
not the case. The ROSAT hardness ratio 1 (HRl) for 
IRXS J175113. 2-293842 is 1.0 ± 0.3 and hardness ratio 
2 (HR2) is 0.54 ± 0.27 (where HR1= (B-A)/(B+A) and 
HR2= (D-C)/(D-|-C), with A=0. 11-0.41 keV, B=0.52- 
2.0 keV, C=0. 5-0.9 keV, and D=0.9-2.0 keV count rate). 
Assuming a power law spectrum with index 2, an inter- 
stellar extinction of 1x10^^ cm~^ (consistent with the 
observed hardness ratio HRl and HR2) and using the 
RASS-measured 0.1-2.4 keV count rate of 0.07 counts 
s~^, we calculate the absorbed 0.1-2.4 keV flux to be ap- 
proximately 1 X 10~^^ erg cm^^ s^^ at the time of the 
ROSAT observation. Assuming Chandra detected zero 
counts and following Gehrels (1986), we take an upper 
limit of 3 counts for a 95 per cent upper limit over the 2 ks 
exposure. This makes the approximate Chandra 0.1-2.4 
keV upper limit for such a source spectrum 5 x 10"^^ erg 
cm~^ s~^, implying a decay in flux of at least a factor 
200. 

For IRXS J174220.8-273736 the detected RASS count 
rate is 0.16 counts s^^, the HRl is 1.00±0.03 and the 
HR2 is 1.00±0.08. For the same assumed source spec- 
trum as above the absorbed 0.1-2.4 keV ROSAT flux 
was 3 X 10^^^ erg cm^^ s^^, implying that the source 
decayed in flux more than a factor of 500. 

Conversely, the bright est sour ce in our survey, 

SAX J1750.8-2900 (NataluccLeTaD 11999), was not de- 
tected in the RASS BSC. This source is a recurrent tran- 
sient LMXB that happened to be in outburst when our 
Chandra observation was obtained but that was appar- 
ently in quiescence during the RASS observations. 

5.3. Sources from pointed ROSAT observations 

We list ROSAT sources discovered in pointed obser- 
vations when they fall inside the covered GBS area. 
Sources found using both the position sensitive propor- 
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TABLE 4 

Sources from pointed ROSAT observations that fall in the 





GBS 


AREA. 






ROSAT name 


Chandra 


a 


S 








(hh:mm:ss) 


(dd:mm:ss) 


IRXH J173547.0-3028511 


# 


43 


17:35:45.53 


-30:29:00.0 


IRXH J173803.6-290706 


# 


31 


17:38:03.50 


-29:07:06.1 


2RXP J173940.7-285116 


# 


115 


17:39:40.80 


-28:51:11.7 


2RXP J174046. 7-283849 


# 


469 


17:40:46.58 


-28:38:50.6 


2RXP J174104.6-2815041 


# 


32 


17:41:04.91 


-28:15:03.4 


2RXP J174133. 7-2840351 


# 


21 


17:41:33.76 


-28:40:33.8 


2RXP J174141.9-2833241 


# 


114 


17:41:41.76 


-28:33:24.2 


2RXP J174249. 9-275028 


# 


785 


17:42:49.95 


-27:50:38.5 


2RXP J174834.7-295730 


# 


230^ 


17:48:35.26 


-29:57:31.9 


2RXP J174834.7-295730 


# 


332 


17:48:35.54 


-29:57:28.8 


2RXP J174928. 4-2919011 


# 


156 


17:49:28.31 


-29:18:59.4 


2RXP J175029. 3-285954 


# 


13 


17:50:29.13 


-29:00:02.3 


2RXP J175041. 2-2916441 


# 


183 


17:50:41.17 


-29:16:44.5 



1 Multiple sources from the ROSAT pointed source catalog are 
consistent with the Chandra position. We list the nearest one. 

2 Those two Chandra sources are 4.8" apart; these are likely two 
detections of the same Chandra source. 

tional counter (PSPC) detector as well as those found 
using the high resolution imager (HRI) observations are 
included (see Table |4]). We exclude pointed observations 
of sources discussed above ( Chandra source numbers 1- 
10). 

Sometimes multiple ROSAT observations of the same 
source region provide slightly different source coordi- 
nates, and thus different ROSAT names, whereas us- 
ing Chandra we find one source. Although it is pos- 
sible that there was indeed more than one source and 
these have faded below the Chandra detection level, 
it is more likely that these multiple ROSAT detec- 
tions are in fact of one and the same source. As 
an example, the six ROSAT sources 2RXP J173826.2- 
290147, 2RXP J173827.4-290138, 2RXP J173827.5- 
290152, 2RXP J173826.7-290205, 2RXP J173827.2- 
290144, 2RXP J173824.0-290146 are probably all re- 
lated to Chandra source #7 (at a=17:38:26.18 and 
29:01:49.4). The nominal angular distance between the 
ROSAT positions and the single Chandra position is 
2.45", 18.4", 16.8", 16.7", 14.9", 28.3", respectively. The 
last ROSAT position was derived from an observation 
that had the source far off-axis (50.8'). Similarly, the 
smallest positional offset was obtained for the observa- 
tion where the source was only 14.8' off-axis. 

6. SUMMARY 

We have started the Chandra Galactic Bulge Survey 
(GBS) with the goals: 

i) to identify quiescent, eclipsing NS and BH X-ray bina- 
ries that are bright enough in the optical or near-infrared 
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to allow phase resolved spectroscopic observations tai- 
lored to measure the compact object mass, 
ii) to constrain binary population synthesis and X-ray 
binary formation and evolution models by means of a 
source number count and a study of the spatial distribu- 
tion of X-ray binaries. 

In this paper we have presented the Chandra source list 
and some properties of the X-ray sources of observations 
covering '^two-thirds (« 8.3 square degrees) of the total 
envisaged survey area of 12 square degrees. The accurate 
Chandra source position will help identify the optical 
counterparts. The 1234 X-ray sources that have been 
discovered so far compares well with the total number of 
w 1650 X-ray sources that we predict we should detect 
in the full 12 square degrees. However, this is of course 
no guarantee that the number of sources per source class 
are close to those we calculated. 

We compared our source list with the source list of 
the RASS. Two BSC RASS sources, IRXS J175113.2- 
293842 and IRXS J174220.8-273736, are not detected in 
our GBS survey, indicating a decrease in flux with factors 
larger than 200 and 500, respectively. Furthermore, we 
compared our Chandra source list with the sources found 
in the catalog of sources derived from pointed HRI and 
PSPC ROSAT observations that fah inside the GBS area 
(see Table 3). 

We also imaged the complete survey area using optical 
observations obtained with the 4m Blanco telescope at 
CTIO. Furthermore, we re-imaged the survey area «30 
times in the r' band to search for (periodic) variable 
sources. Results of these campaigns will be presented 
in forthcoming papers. 

7. OUTLOOK 

We are in the process of obtaining optical spectro- 
scopic and photometric observations of all optical pro- 
posed counterparts in the Chandra error circle to clas- 
sify the X-ray sources. This step is crucial to achieve our 
science goals outlined above. 

A full discussion of all known information on all the X- 
ray sources is beyond the scope of this paper and we defer 
the discussion of the other fainter Chandra sources to a 
forthcoming paper, where they will be discussed together 
with optical photometric and spectroscopic information. 
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